ON THE GREEN'S MATRICES OF STRONGLY PARABOLIC 
SYSTEMS OF SECOND ORDER 



SUNGWON CHO, HONGJIE DONG, AND SEICK KIM 

Abstract. We establish existence and various estimates of fundamental ma- 
trices and Green's matrices for divergence form, second order strongly par- 
abolic systems in arbitrary cylindrical domains under the assumption that 
solutions of the systems satisfy an interior Holder continuity estimate. We 
present a unified approach valid for both the scalar and the vectorial cases. 



1. Introduction 

In this paper, we study Green's matrices (or Green's functions) of second order, 
strongly parabolic systems of divergence form 

N N n 

(P) A^^'-^ ^ i?.„(v4r/(i,2:)i?,,zi^"), z = l,...,iV, 

j = l j = l a,/3=l 

in a cylindrical domain U = M. x fl where is an open connected set in E" with 
n > 1. By a Green's matrix for the system jP]) inU = M x we mean an x iV 
matrix valued function Tij{t, x, s,y) {x,y € Q and t,s G M) which satisfies the 
following (see Theorem 12. 71 for more precise statement): 

N 

lim Tij{t,x,s,-) ^SijSxi-), 

t — >s-\- 

Py(.,.,s,y) = on dU^Rxdn, 

where dij is the Kronecker delta symbol, i5(s,j^)(', •) and Sx{-) are Dirac delta func- 
tions. In particular, when U = M"+^, the Green's matrix (or Green's function) is 
called the fundamental matrix (or fundamental solution). 

We prove that if weak solutions of jP]) satisfy an interior Holder continuity 
estimate (see Section [^751 for the precise formulation), then the Green's matrix exists 
in l/( and satisfies several natural growth properties; see Theorem 12.71 Moreover, 
when U — we show that the fundamental matrix satisfies the semi-group 

property as well as the upper Gaussian bound of Aronson [T] ; see Theorem 12.111 
The method we use does not involve Harnack type inequalities or the maximum 
principle, and works for both the scalar and the vectorial situation. Moreover, we 
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do not require the base fl of the cyhnder U = R x to he bounded or to have a 
regular boundary. In the scalar case (i.e., N = 1), such an interior Holder continuity 
estimate for weak solutions is a direct consequence of the celebrated theorem of Nash 
[24] (see also Moser ,23 ), and thus, we in particular prove the existence and growth 
properties of Green's functions for the uniformly parabolic equations of divergence 
form with bounded measurable coefficients in arbitrary cylindrical domains; see 
CoroUarv 12.81 If n = 2 (i.e., fl C M^) and the coefficients A"^f^{t,x) of the system 
(|P|) are independent of t, then we will find that the Green's matrix exists in any 
cylindrical domain and satisfies the growth properties stated in Theorem \2.7\ see 
CoroUarv 12.91 Also, if > 1 and the coefficients Afj^{t,x) of the system dP]) are 
uniformly continuous in x and measurable in t or belong to the class of VMO, then 
we have the same conclusion; see Corollary 12.101 

Let us briefly review the history of work in this area. Fundamental solutions of 
parabolic equations of divergence form with bounded measurable coefficients have 
been a subject of research for many years. The first significant step in this direc- 
tion was made by Nash [24| who established certain estimates of the fundamental 
solutions in proving local Holder continuity of weak solutions. Aronson [T] proved 
Gaussian upper and lower bounds for the fundamental solutions by using the para- 
bolic Harnack inequality of Moser [23] . Fabes and Strook ^ showed that the idea 
of Nash could be used to establish Aronson's Gaussian bounds, which consequently 
gave a new proof of Moser's parabolic Harnack inequality. In [2], Auscher gave a 
new proof of Aronson's Gaussian upper bound for the fundamental solution of para- 
bolic equations with time independent coefficients, which carries over to the case of 
a complex perturbation of real coefficients. Recently, it is noted in [14] that Aron- 
son's upper bound is equivalent to the local boundedness property of weak solutions 
of strongly parabolic systems. Green's functions of elliptic equations of divergence 
form in bounded domains have been extensively studied by Liftman, Stampacchia, 
and Weinberger [50j and Griiter and Widman , whereas the Green's matrices of 
the elliptic systems with continuous coefficients in bounded domains have been 
discussed by Fuchs [TT] and Dolzmann and Miiller [S]. Very recently, Hofmann 
and Kim |13j gave a unified approach in studying Green's functions/matrices in 
arbitrary domains valid for both scalar equations and systems of elliptic type by 
considering a class of operators L such that weak solutions of Lu = satisfy an 
interior Holder estimate. Some parts of the present article may be considered as a 
natural follow-up of their work in the parabolic setting. Readers interested in the 
construction of fundamental matrices for parabolic systems in nondivergence form 
with Holder continuous coefficients are asked to refer to Eidel'man Friedman 
[To] , or Ladyzenskaja et al. Also, we would like to bring attention to a paper 
by Escauriaza 7] on the fundamental solutions of elliptic and parabolic equations 
in nondivergence form. 

The organization of this paper is as follows. In Section |21 we define the property 
(PH) for parabolic systems in terms of Morrey-Campanato type Holder estimates 
for weak solutions. As we mentioned earlier, in the scalar case, such a property is 
a direct consequence of the interior Holder continuity estimates by Nash [J^ and 
Moser [23l. We also present some other examples of parabolic systems satisfying 
the property (PH), including the almost diagonal systems and the systems with co- 
efficients from the class VMO^ (see below for the definitions) . We close the section 
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by presenting our main theorems. In Theorem 12.71 we state the existence, unique- 
ness, and properties of Green's matrices in an arbitrary cyhndrical domain for the 
parabohc systems satisfying the property (PH) . Theorem 12.111 estabhshes Gauss- 
ian bounds and semi-group properties for the fundamental matrices of parabohc 
systems satisfying the property (PH) . In Section [3l we prove Theorem 12.71 in the 
entire space assuming that the property (PH) holds globally. In Section HI we give 
the proof for Theorem 12.71 in the general setting. Finally, we prove Theorem 12.111 
in Section [S] 



2. Preliminaries and main results 

2.1. Basic Notations. We use X = {t,x) to denote a point in M"+i = R x R" 
with n> \] X = (xi, . . . , Xn) will always be a point in M". We also write Y ~ (s, j/), 
-'''o = [to^xo), etc. We define the parabolic distance between the points X = (i,a:) 
and Y = {s,y) in M"+i as 

\X - Y\p := max(v/|i - s|, \x - y\), 

where | • | denotes the usual Euclidean norm. We usually use lA to denote an open 
set in R"+^ and Q. to denote an open set in R". We define 

Br{x)^{y^W' ■.\y^x\<r} 

and use the following notations for basic cylinders in 

Q-{X) = it-r^,t) X Brix), 
QtiX) = {t,t + r^) X Br{x), 
Q^{X) = (t~r^,t + r^) X Br{x). 

Note that Qr{X) = {Y e W'+^ : \Y - X\p <r}. We use dpQ-{X) and dpQ+{X) 
to denote the parabolic forward and backward boundaries of Q~{X) and Q^{X), 
respectively; i.e., 

dpQ-{X) = {t- r^,t) X dBr{x) U {t - r^} x Br{x), 
dpQ+{X) = (t, i + r^) X dBr{x) U {t + r^} x Br{x), 

where dBr{x) and B^ix) denote the usual boundary and closure of Br(x) in M"; 
i.e., dBrix) = {yeW:\y-x\^r} and Br{x) = {y G K" : |y - a;| < r}. 

For a given function u = u{X) — u{t,x), we use DiU for du/dxi while we use 
ut (or sometimes Dtu) for du/dt. We also write Du (or sometimes D^u) for the 
vector (Diu, . . . , Dnu). If u is a function in a set Q C M"+^, we denote 

\u{X)-u{Y)\ , , , 

[mJcm,m/2(q) := sup — — — — , where ^ e (0,1J. 
x^Y \X -Y\p 
X,YeQ 

For a Lebesgue measurable set Q C (resp. S C M"), we write \Q\ (resp. \S\) 

for the Lebesgue measure of the set Q (resp. S). For u G L^{Q) (resp. u £ L^{S)) 
we use the notation j-Q "u- = Jq u (resp. j-gU = -rir Jgu). 
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2.2. Function Spaces. We follow the notation of 20J with a slight variation. For 
U C M"+^, we write U{ti)) for the set of all points (to, x) in U and liU) for the set 
of all t such that U (t) is nonempty. We denote 

\Mu = \\Du\\l2(u^ +CSSSUp|lu(i,-)|1^2(;^(t)). 

tei{U) 

For u — {u^,. . . , u^), we write |||m|||^ := J2^=i 

\Ml = +eSSSUp||M(t,-)||i2(;^(t)). 

te/(W) 

In the rest of this subsection, we shall denote by Q the cylinder (a, b) x il, where 
— oo<a<6<oo and is an open connected (possibly unbounded) set in R". We 
denote by W^^iQ) the Hilbert space with the inner product 



and by (Q) the Hilbert space with the inner product 

(^i,^;)^yl,l(Q^ uv + DkuDkV + / UtVt- 

We define V2{Q) as the Banach space consisting of all elements of W2'^{Q) having 
a finite norm ||u||v'2(q) ■= III'^IIIq ^^'^ V^z'^iQ) the Banach space consisting of all 
elements of VziQ) that arc continuous in t in the norm of L^(f7), with the norm 



II"IIIq = (^ll^w|li2(Q) + maxJ|u(i,-)||i2(o)) 



1/2 



The continuity in t of a function u{t, x) in the norm of L'^{Q) means that 
lim \\u{t + h,-)- u{t, •)|| = 0. 

The space V2^'°(Q) is obtained by completing the set W2'^{Q) in the norm of VziQ). 

We write C^p{Q) for the set of all functions u e C°°{Q) with compact supports 
in [a,b] x ft while C^iU) denotes the set of all infinitely differentiable functions 
with compact supports in U. We denote by Wl'^{Q) and Wl''^{Q) the closure 
of C^((5) in the Hilbert spaces Wl'^{Q) and W2^{Q), respectively. We define 
V2{Q) ■= V2{Q)nW2-°{Q) and F2^'°(Q) = 1^2^'°(Q)nM^2 '°(Q). It is routine to check 
that VziQ) and V2 '^{Q) are subspaces of the Banach spaces V2{Q) and V2^'°(Q), 
respectively. By a well known embedding theorem (see e.g., [20', §11.3]), we have 

(2.1) lklL2+4/„(Q) < C(n)|||w|||Q ^ueUQ)- 

When is an infinite cylinder (i.e., (a,oo) x ft, (—00,6) x il, or (—00,00) x fl), we 
say that u € V2{U) if m e VziUr) for all T > 0, where Ut ■.= U^ (-T, T) x M", and 
< 00. Similarly, we say that u e V2{U) (resp. ^^'"(W), V'2^'°(Z^)) if u e V2(ZiT) 
(resp. V'2^'°(Wt), V"2^'°(Z^t)) for aU T > and < 00. Also, we write u € C^piU) 
if w € C^p{Ut) for all T > 0. For a sequence {ufcl^^ in V2{U), we say that Uk ^ u 
"very weakly" in VziJA) for some u S V2(Z^) u weakly in W2'^{Ut) for all 

T > 0. If {ufc}^^ is a bounded sequence in V2(Z^) (resp. VziUj), then there exists 



GREEN'S MATRICES 5 

a subsequence {ukj}fLi ^ {"felfe^i ^^'^ ^ ^ V2{U) (resp. V2(W)) such that Ukj u 
"very weakly" in V2{U) (see Appendix for the proof). 

The space W^''^{U) {1 < q < oo) denotes the Banach space consisting of functions 
u € L''{U) with weak derivatives DaU £ L'^iU) (a = 1, . . . , n) with the norm 

We write u S L^iU) if u S L°°{U) has a compact support in U. 

2.3. Strongly parabolic systems. Throughout this article, the summation con- 
vention over repeated indices are assumed. Let C — dt ~ L he a. second or- 
der parabolic operator of divergence type acting on vector valued functions u = 
{u^, . . . , u^Y" ^ 1) defined on an open subset of in the following way: 

(2.2) Cu^ut- Lu:=ut~ Da{A°''^ Dpu), 

where A"'^ (a,/3 = 1,... , n) are N hy N matrix valued functions with components 
(A^^)f^-^j^ defined on M"+^ satisfying the strong parabolicity condition, i.e., there 
is a number A > such that 

N n 

(2.3) Af{X)epC.>\\i\^:=xY,Y.\^l.\'' 

i=l a = l 

We also assume that A'^^ are bounded, i.e., there is a number A > such that 

N n 

(2.4) ^ |Ar/(A)|2<A2 VAeM"+i. 

i j' — 1 a,/3— 1 

If we write (j2.2[) component-wise, then we have 

{Cuf ^u\- D^{AfDpu^) - 1, . . . , TV. 
The transpose operator *£ of C is defined by 

where = (yl^")^ (i.e., *Aff = A^")- Note that the coefficients *Aff also 

satisfy (|2.3p and (|2.4p with the same constants A, A. 

2.4. Weak solutions. For f,g^ G L\^^{U)'^ (a = l,...,n), we say that u is a 
weak solution of Cu = f + DaQ^ U if u € V2{U)^ and satisfies 

(2.5) - / v}^ + I AfDpuW^q^^= [ fcj,^- [ g\D^^^ yct>eC^iUf. 

Ju Ju Ju JU 

Similarly, we say that m is a weak solution of *£tt — f + DaOa in if m G V2(^)^ 
and satisfies 

(2.6) / u^^l + [ 'AfDpuW^q^^= I rcj>^- I g\D^^^ y4>£C^iU) 
Ju Ju Ju Ju 

By a weak solution in V2'^'"((a, b) x Vl)^ (— oo < a < b < oo) of the problem 



N 



(2.7) 



Cu = f 
u{a, •) = g, 
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we mean a function u{t, x) that belongs to V2 ' ((a, b) x O) and satisfying for all 
ti in [a, 5] the identity 



i1 J a J fl 



Hi 



f^' + / gV^(a, •) V0 e b) x 17)^ 



Similarly, by a weak solution in V2 ' (('ij ^) ^ of the (backward) problem 

uib, •) = g, 



(2.9) 



we mean a function u{t, x) that belongs to ^2 '"(('^1 ^) ^ ^^^^ satisfies for all ti 
in [a, 6] the identity 

MXii,-)0'(ii,-)+ / / [ 'A^fDpuW^q^' 

b 



(2.10) = / / / .gVX^-) V0eC?°p((a,5)xn) 

We say that m is a weak solution in V^ '^ {\a^co) x f2)^ of the problem (|2.7p if 
M is a weak solution in V'2^'°((a, 6) x £7) of the problem ()2.7|) for all & > a and 
lll'"lll(a,oo)xo < oo- Similarly, we say that m is a weak solution in y2^'*'((— cx), 6) x f2)^ 
of the problem (|2.9p if it is a weak solution in V'2^'°((a, 6) x il) of the problem (|2.9p 
for all a < 6 and |||M|||(-oo,b)xo < oo- 

Lemma 2.1. Assume that H. is an open connected set in M", / G iJ!°(]R x il)^ , 
and g G i^(f2)^. Then, there exists a unique weak solution in V'2^''^((a, cxd) x ft)'^ 
(resp. l/2^'°((-oo,6) x fl)'^ ) of the problem fSJl) (resp. !^ ). 

Proof. See e.g., [20, §111.4] for the existence and PO', §111.3] for the uniqueness. We 
point out that the proof does not require the boundedness of fl. □ 

2.5. Property (PH). We say that the operator C satisfies the property (PH) if 
there exist fiQ G (0, 1], Rc E (0, 00], and Co > such that all weak solutions u of 
Cu = in ~ Q]^{Xq) with R < Rc satisfy 

(2.11) / \Duf < Co ( -Y^^''° [ \Duf \/0<p<r<R. 

Similarly, we say that the operator *£ satisfies the property (PH) if all weak solutions 
u of *Cu = in Qjij Q^{Xo) with R < Rc satisfy 

(2.12) / \Duf < Co ( ^Y^^''° [ \Du\^ yO<p<r<R. 
jQt ^""^ JQt 

Now, we present examples of strongly parabolic systems satisfying the property 
(PH). Lemma l2.2l states that almost diagonal parabolic systems satisfy the property 
(PH) with Rc = 00 and Lemma [2.31 shows that parabolic systems with coefficients 
which belong to VMO^; (see below for the definition) satisfy the property (PH) for 
some Rc < 00. In [TBI, it is shown that if n = 2 and the coefficients of £ are 
independent of t, then £ and *£ satisfy the property (PH) with Rc — 00. 
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Lemma 2.2. Let {o,"^ {X))^ be coefficients satisfying the following conditions: 
There are constants Aq, Aq > such that for all X G M."''^^ 

n 

a"''WC/3ea > Ao 1^1' VeeM"; ^ \a"^{xf<Al. 

a, (3=1 

Then, there exists Sq = EqIu, Xq, Aq) > such that if 

N n 

^'iX)--= E E \A?f{X)-a-^{X)S., 

= l a, (3=1 

where Sij is the Kronecker delta symbol, then the operator C associated with the coef- 
ficients A"j^ and its transpose */! satisfy the property (PH) with /io = /xo(n, Aq, Aq), 
Gq — Co{n, N, Ao, Ao), and Rc — oo. 

Proof. See e.g., P3, Proposition 2.1]. □ 
For a measurable function / = f{X) = f{t,x) defined on M"+^, we set 

ujsif) sup sup ] [ I \f(jj,s) - fx,r{s)\ dyds V(5 > 0, 

Jf = (t,a;)eR"+i r<(5 \Ulr\^)\ Jt-r^ J B,.{x) 

where fx,r{s) = j-B,{x) ')• ^^^^ / belongs to VMO^; if lims^oCLisif) = 

0. Note that VMO^; is a strictly larger class than the classical VMO space. In 
particular, VMO^^ contains all functions uniformly continuous in x and measurable 
in t; see [TO] . 

Lemma 2.3. Let the coefficients of the operator C in (12. 2p satisfy the conditions 
(|2.3p and (12. 4p . // the coefficients belong to VMO^;, then the operators C and 
satisfy the property (PH). 

Proof. Let u, be a weak solution of Cu — in Q^{Xo) with R < Rc, where Rc 
is to be chosen later. First, note that as a weak solution of a linear strongly 
parabolic system, u satisfies the following improved integrability estimates for some 
p = p{n, A, A) > 2 (see [HI Theorem 2.1]): 

(2.13) (i \DuA <cS \Du\' Vr<i?/2. 



\jQr(Xo) J ''Q2riX0) 

For a given r < R/2, denote A^f = A-^l^^^{t) = j^BA^o) ^"'^(*' ■)■ We decompose 
u = V + w, where w is the weak solution in V2''^{Q~ (Xq))^ satisfying 
(2.14) wt - Dc.{Af^^^Dpw) = Dc.{{A'^P - Af^^^)Dpu) in Q- {X^) 

with zero boundary condition on dpQ^{XQ). By using w itself as a test function 
in (|2.14p and then using Holder's inequality, (|2.4p . and p.l3p . we derive 

\ 2/9 / „ \ 2/p 



/ \Dw\^<c( f \A-A,,,r\'] { I iDul" 



(2.15) < CA2(?-i)/«a;^(A)2/'? I \Duf , where q = 2p/{p - 2). 

JQ2M0) 

Observe that v — u — w is a. weak solution of 

vt - Da, {A'^^^,{t)Dpv) = in Q;{Xo). 
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Since is a weak solution of a strongly parabolic system with coefficients indepen- 
dent of the spatial variables, v possesses the following interior estimate: 

(2.16) / \Dvf<c(-Y^^ [ \Dvf yp<r, 

for some constants C — C{n, N, X, A) (see Appendix). Then, by combining (|2.15p 
and (|2.16p . we see that for all p < r < R/2, we have 

\Du\^ < C (^y^^ f iDul^ + CiOriAf/" f 

If we choose Rc sufficiently small so that Cujr^ (A)^/' is small, then by a well-known 
iteration argument we obtain (|2.1ip (see e.g., [T51 Lemma 2.1, p. 86]). Therefore, 
we have proved that £ satisfies (PH) . The proof that *£ also satisfies (PH) is similar 
and left to the reader. □ 

2.6. Some preliminary lemmas. 

Lemma 2.4. There exists a constant C = C{n, A) > such that if u is a weak 
solution of Cu — f in Q~j^ — Q^{Xq), then 



where Un — j^^- u. Similarly, if u is a weak solution of ^Cu = f in Q^, then 
where Ur = -f„+ u. 

Proof. See e.g., |271 Lemma 3]. □ 

Lemma 2.5. Let u G _L'^((52ij(^o))^ '^nd suppose that there exist positive constants 
fj, G (0, 1] and Nq such that 

I <iVo'''"+'+'^ VXgQ^(Xo) VrG(0,i?), 

"'Qr (X) 

where ux,r — fq- (^x) ^^s^' ''■^ Holder continuous in Q^(Xq) and 

[m]cm,^./2(q^(Xo)) < C{n,N,^i)No. 
The same is true with in place of everywhere. 

Proof See e.g., pH Lemma 4.3]. □ 

Lemma 2.6. Assume that the operator C satisfies the property (PH). Then, all 
weak solutions u of Cu = in = Q^(Xq) with R < Rc satisfy 

^1/2 

(2-17) IIHL^(Q-^jXo))<c' ' 



\jQR{Xa) 



where C = C{n, N, X, A, fiQ, C'o) > 0. Moreover, all weak solutions u of Cu = in 
Qr = Qi?,(^o) with R < Rc satisfy 

(2.18) ML^iQ,)^ ^R_r)U2)/p Ml^Qr) VrG(0,i?) Vp > 0, 
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where Cp — Cp{n, N, A, A, /io, Co) > 0. ^ similar statement is also true for *£ with 
replaced by everywhere. 

Proof. Let M be a weak solution of Cu = in Q^{Xo), where R < Rc. For any 
X e Q~^^^{Xq) and r < i?/4, it follows from Lemma [2.41 property (PH), and the 
energy inequality (see e.g., [20l §111.2]) that 

/ \u - ux.rf < Cr^ f \Du\'^<Cr'^{r/R)''+'^''°f \Duf 
<C{r/R)''+'^+^f'« [ \uf 



(2.19) < C'j.»+2+2mo^-2mo i 

Jq 



\u\\ 



Note that (j2.19|) also holds for any X € Q^^yjl^o) and any r S [i?/4, R/2), because 



/ 



iQ7(X) JQ7(X) JQr(Xo) 

Therefore, by Lemma |2.5[ we obtain 



(2.20) Ma™,.o/.(Q,^,(.,)) < ^^-^^°/ 



Then, p.l7p is an easy consequence of p.20p and a well known averaging argument 
(see e.g., [H]). We point out that it actually follows from (|2.19p 

(2.21) ,^,<Ci?-^-/ |Z?«r. 

For the proof that ((2Tf|) implies ((2l8)) . we refer to [HI pp. 80-82]. □ 
2.7. Main results. We now state our main theorems. 

Theorem 2.7. Let = M x f2, where is an open connected set in R". Denote 
dx = disi{X,dU) = <l\si{x,dn) for X = {t,x) € U; we set dx = oo if Vl = R". 
Assume that operators C and *£ satisfy the property (PH). Then, there exists a 
unique Green's matrix T{X, Y) = r(t, x, s,y) on hi hi which is continuous in 
{{X, Y) hi X hi : X ^ Y}, satisfies T{t, x, s,y) = for t < s, and has the property 
that T{X, .) is locally integrable in hi for all X £ hi and that for all f € C^{hi)^ , 
the function u given by 

(2.22) u{X) := / T{X,Y)fiY)dY 

Ju 

belongs to V2 '^{hi)^ and satisfies Cu = f in the sense of (j2.5p . Moreover, T 
satisfies 



(2.23) [-T,k{;Y)^l + AfDpr,ki;Y)D^(j)^) = ^\Y) V0 G CT{U) 

and for all rj G C^'(hi) satisfying rj = I on QriY) for some r < dy, we have 

(2.24) {l^^)Ti;Y)eV,'''{U) 



N 
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Furthermore, for all g G L'^{il)^ , the function u{t,x) given by 

(2.25) u{t,x):= / T{t,x,s,y)g{y)dy Wx e n yt > s 

Jn 

is the unique weak solution in V2''^{{s,oo) x fl)'^ of the Cauchy problem 

(2.26) 1^7 = 

and if g is continuous at xq € in addition, then 

(2.27) lim [ r{t,x,s,y)g{y)dy^g{xo). 

Denote dx ■= miii(dx, Rc) for X Cz U. Then, T satisfies the following estimates: 

(2.28) l|r(-,r)IL.+4/„(^,\Q^(^)) + |||r(.,r)||^Q^(^) < Cr-"/^ Vr < dy, 

(2.29) \\T{X, ■)h.+.,.^u\QAX)) + ll|r(^, ■)\lu\QAX) < Gr^'^ Vr < dx, 

(2.30) ||r(.,y)|Up(Q^(y)) <Cpr-"+("+2)/f Vr<Jy Vpe[l,ii±^), 

(2.31) ||r(X,.)||LP(Q„(x)) <Cpr-"+("+2)/f Vr<dx Vp £ [1, i^), 

(2.32) |{X e W : |r(X,y)| > t}| < Ct-("+2)/" Vr>(dy/2)-" VF e Z^, 

(2.33) |{y e : |r(X, r)| > r}| < Cr-("+2)/" Vt > (rfjc /2)-" VX G Z^, 

(2.34) ||Z)r(.,r)|U.(Q^(y)) <Cpr-"-i+("+2)/f Vr<dV Vp G [1, ^), 

(2.35) ||OT(X,.)||LnQ.m) <^P^"""'+'"+''^'' Vr<dx Vp G [1, ^), 

(2.36) \{X(^U: \DJ^(X, Y)\ > t}\ < Ct^"^ Vt > (Jf/2)-" VF G Z^, 

(2.37) \{Y eU : \Dyr(X, Y)\ > t}\ < Ct-"^ Vt > (djc /2)-" VX G Z^, 

(2.38) |r(x,r)| <qx-y|;" o < |x-y|p < iniax(Jx,dV), 

w/iere C = C(n, TV, A, A,/io,Co) > and Cp = Cp{n, N, \, A, fio,Co,p) > 0. 

Corollary 2.8. If N ~ \ or if N > 1 and the coefficients of C satisfies the as- 
sumption of Lemma \2.'A then there exists a unique Green's function in a cylindrical 
domain satisfying the properties in Theorem \2. 7| with Rc — oo. 



Proof. See Lemma [2T2l □ 

Corollary 2.9. If n = 2 and the coefficients of C are independent oft, then there 
exists a unique Green's matrix in a cylindrical domain satisfying the properties in 
Theorem \2. 7| with Rc ^ oo. 

Proof. See P^l for the proof that C and *£ satisfy the property (PH). □ 

Corollary 2.10. If N > 1 and the coefficients of £ belong to VMO^r (e.g. the 
coefficients are uniformly continuous in x and measurable in t or belong to VMO ), 
then there exists a unique Green's matrix in a cylindrical domain satisfying the 
properties in Theorem\2.7\ with Rc > 0. 



Proof. See Lemma [^751 □ 
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Theorem 2.11. Assume that operators C and satisfy the property (PH) and let 
T(t, X, s, y) be the fundamental matrix for C inlA = W^^^ as given in Theorem \2.7\ 
Then, we have for all t > s and x,y Cz M" 

(2.39) \Tit, X, s, y)l^ < C{t - s^^/' exp{-K\x - y\^/{t ~ s) + ^{t - s)/RI}, 

where C = C{n, N, X, A, ^o,Co), k = k(A, A),7 — 7(71, A^, A, A, /io, Cq) > 0, and 
I ' lop denotes the operator norm. In particular, when Rc — 00, we have the following 
usual Gaussian bound for all t > s and x,y £ K"; 

(2.40) \r{t, X, s, y)lp < C{t - s)-"/2 exp{-«|x - y\y{t - s)}. 
Moreover, for t > s the following identities hold: 

(2.41) r{t,x,s,y)= r{t,x,r,z)T{r,z,s,y)dz Va;,?/eM" Vre(s,i), 

(2.42) / T{t,x,s,y)dy = I Vx € R", 
where I is the N by N identity matrix. 

Remark 2.12. In fact, ((2:4T|l also holds for n^W"-; see Section [521 



3. Proof of Theorem 12.71 when il = M" and Rc = 00 

In this section we shall prove Theorem l2 . 7l under the assumption that = R" and 
that the operators C and *£ satisfy the property (PH) with Rc — 00. Consequently, 
we have U = R"+^ and dx = dx — 00 for all X £ R"+^. Throughout this section, 
we employ the letter C to denote a constant depending on n, N, A, A, /xq, Cq, and 
sometimes on an exponent p characterizing Lebesgue classes. 

3.1. Averaged fundamental matrix. Our approach here is an adaptation of 
that in Hofmann-Kim [13] , which in turn is partly based on the method by Griiter 
and Widman [H]. Let F = (s, y) S R"+i and 1 < A: < TV be fixed. For each p > 0, 
fix So G (~c», s — p^). We consider the problem 

(3.1) I ^" = ^V(n^'= 

[ m(so,-) = 0, 

where is the k-th unit vector. By Lemma [2.11 we find that the problem p.ip 
has a unique weak solution Vp = Vp.y^k in 1^2^'°((so, 00) x R")^. Moreover, by the 
uniqueness, we find that Vp does not depend on the particular choice of sq and we 
may extend Vp to the entire R"+^ by setting 

(3.2) Vp = on (-oo,s- p2) X R". 
Then Vp £ V^'° and satisfies for all ti > s the identity 



(3.3) / vl<j>\t,,-)^ r [ vlcj>\+ r f AfDpv^Dc^d^^^i 

JR" J-ooJR" J-ooJR"- JQ] 

for all 4> £ C^(R"+^)^. We define the averaged fundamental matrix r''(-,y) 
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Next, for each / e if (M"+i)^, let us fix to such that / = on [to,oo) x 
We consider the backward problem 



(3.4) 



'Cu = f 

M(t0,-) =0. 



Again, by Lemma [2A] we obtain a unique weak solution u in V^^'°((— oo, to) x M."')^ 
of the problem (|3.4p and we may extend u to the entire by setting m = on 

(io,oo) X K". Then, u € F2^'°(M"+i)^ and satisfies for all ti the identity 



(3.5) / uVXii,-)+ / / 



for all ^ e C,^p(M"+i) 



Lemma 3.1. For Vp and u constructed above, we have 

(3.6) |K|||r„+. <c|Q;(r)p^, 

(3.7) IHk..i<C 11/11^^^^^^^^, 



(3.8) / ■ / = / 

JR"+1 JQ 



Proof. The energy inequality (see e.g., [20l §111.2]) together with (|2.ip yields (|3.6p 
and p.7p . The identity p.Sp follows from p.3p and p.Sp accompanied by standard 
approximation techniques (see e.g., §111.2] or f5P, §VI.l] for details). □ 

3.2. L°° estimate for averaged fundamental matrix. Let u g ¥2'^^ iW^^^)^ be 
constructed as above with / e Lf{W^+^)^. Fix Xo e W+\ R>0,Xg Qr{^o), 
and r g (0, i?]. We decompose u = ui + U2, where U2 is the unique weak solution 
in V^''^ {Qt {X)) of *£it2 = / in Q^{X) with zero boundary condition on dpQ^{X). 
Then, Ui = u — U2 satisfies *£iti = in Q^{X), and thus, for < (5 < r, 

|Dtt|' < 2 / |i?i6i|' + 2 / \DU2\^ 

-(X) JQt(x) jQtm 



< C(<5A)"+"'" / \Duif + 2[ \Du2f 

JQtiX) JQtiX) 



(3.9) <CiS/rr+^^° \Duf + C \Du2\\ 

Jq+(x) JQtiX) 

For a given p> {n + 2)/2, choose pg ^ ((^^ + 2)/2,p) such that 

^1 := 2 - (n + 2)/po < Mo. 

As in (13.71). we have 



(3-10) I \DU2\' <C\\f\\\_^ < Cr-+'^^ ll/llLo(Q+(x)) 

Combining (|3J)l with ((3?T0)) . we get for all 5 < r < i?, 

/ \Du\' < C(5/r)"+2^« / |i?«|VCr"+2A'^ II/IILo(q+(x)) • 
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Then, by a well known iteration argument (see e.g., |15l Lemma 2.1, p. 86]), 



(3.11) / \Du\^ <C{r/RT+^^'^ ( |i?M|' + Cr"+2^H|/|| 
By Lemma [231 p.lip . and Holder's inequality, we get 



2 

LPO(Q+(X)) ' 



^^^^^ l« - ^r\' < WDuWl.^Q^^^^^^ + ||/|iL(Q+(X))) 



Then, from Lemma 12.51 and (|3.7p , it follows that 

Mc^'l'^'l/2(Q+(Jfo)) - " ^^MI^'"IIl2(r„ + 1) + ||/||iPO(R'«+i) 



(3.12) 

By Holder's inequality, ^J^, and (|377l) . 

(3.13) < Ci?2 ll/f ,„+4 

By p.l2p . ()3.13p . and a standard averaging method (see e.g., [Il|), we obtain 



"lli=°('3i./2(-^o)) - ^"Jc''i.^'i/2(Q+(Xo)) II"I!l2(Q+(Xo)) 

2 

LP 



(3.14) < Ci?-" 11/11^^^^^^^^ + CR'^^ ||/|lL„(K„.i 

In the remaining part of this subsection, we shall assume that / is supported in 
Q~l^{Xo). Recall that p> {n + 2)/2. We apply Holder's inequality in (|3.14|) to get 

(3-15) II«IIl==(q+,,(Xo)) < II/IIl.(qS(Xo)) • 

If QpiY) C Q+/2(Xo), then dH^HD together with yields 

Jq+(Xo) JQpiY) 

for any p > {n + 2)/2. By duahty, it follows that if Qp(V) C Q^^^{Xo), then 

(3-16) II'^pIIl.(q+(Xo)) < Vq G [1, (n + 2)/n). 

Consequently, we obtain the following pointwise estimate for the averaged funda- 
mental matrices, which is the main result of this subsection. 

Lemma 3.2. Let X = (t,x), Y — (s,y), and assume X ^Y. Then 

(3.17) \T''iX,Y)\<C\X-Y\;" \fp<\X~Y\^/3. 

Proof. Denote d = \X - Y\^ and let Xq ^ {s - Ad'^,y), r = d/3, and R = 20r. It is 
easy to see that (recall p < r) 

q;{y) c Q+/2(Xo), q;{x) c q+(Xo). 

Moreover, by p.3p . Vp — Vp.y.k is a weak solution of £m = in Q^{X). Therefore, 
by Lemma [^751 and p.l6p . we have 

\vp{X)\ < Cr—^ II«pILhq.-(x)) < Cr-'^'" II'^pIIlhq+(Xo)) < 
The lemma is proved. □ 
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3.3. Construction of the fundamental matrix. In next two lemmas, we derive 
LP estimates uniform in p > for averaged fundamental matrices V^^-^Y) and their 
spatial derivatives DT''{-, Y). 

Lemma 3.3. For any Y E and any p > 0, we have 

(3.18) I {X e M"+i : \D^TP{X, Y)\ > t} \ < Ct^"^ Vt > 0, 

i.e., |£'r''(-,y)| has a zweaA;-L^"+^^/("+-'-' estimate in M"+-'-. Moreover, for any 
R > 0, the following uniform estimate for DT'^{-, Y) holds: 

(3.19) ||OT^(.,y)||^,(Q^(^))<C,i?— Vp>0 Vpe[l,^). 
Proof. We claim that for all i? > 0, we have 

(3.20) \lvX^.r\Q,iY)<CR-'' Vp>0. 

To prove the claim (|3.20p . we only need to consider the case R > 6p. Indeed, if 
R < 6p, then ([32]) yields 

i^pi|..+i\Q,(y) < i^piM"+i < Cp-'^ < CR-". 

Fix a cut-off function ( e C^{Qr{Y)) such that 

C=lonQ^/2(y), 0<C<1, \DC\<CR-\ |Ct| < Ci?-^. 

By following a standard proof of the energy inequality (see e.g., [20l §111.2]), we 
derive from ((3TT|l and ([STT]) that 

sup / {l-Cr\v,it,-f+ [ il-0'\Dv/ 

<C [ (\DC\' + \{l-OC,\)\v/ 

(3.21) <CR-^ [ \X -Y\~^"dX <CR-". 

J{B./2<\X-Y\p<R} 

So, we have proved the claim (|3.20p . Note that we also obtain from (|3.2ip 

(3.22) |||(1-C)^pIIIr™+i <Ci?-" Vp>0. 

Now, let Ar = {X e R"+i : \DvP{X)\ > t} and choose R ^ T~^/i^+^). Then, 

\Ar\QRiY)\<T-^ [ iDvPf^CT-"^ 
JAr\QR{Y) 

and \ArnQR{Y)\ < CR'^+^ = Ct'"^ . We have thus proved To show 

(I3.19p . we first note that, for any t > 0, we have 

JQaiY) JQR{Y)n{\Dv''\<r} J QR{Y)n{\DvP\>T} 

(3.23) <tP\Qr\+ I \DvX- 

J {\DvP\>t} 
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By using (|3.18p and the assumption p £ [1, ^^), we estimate 

f \Dvp\'' = f ptP-^ \{\DvP\ > max(<,r)}| dt 

J{\Dvp\>t} Jo 

/"OO 



(3.24) ^c{l-p/{p~^)) rP^("+2)/(«+i) . 
By setting r = in ([3?23l) and ((3?24|) . we get 

from which (j3.19p foUows. The lemma is proved. □ 

Lemma 3.4. For any Y E and any p > 0, we have 

(3.25) |{X e M"+i : |r''(X,y)| > r}| < CT-("+2)/n Vr > 0, 

i.e., |r''(-,y)| has a wea/c-L^"^^-'/" estimate. Moreover, for any R> 0, the follow- 
ing uniform estimate for T'^{-, Y) holds: 

(3.26) ||r''(,y)||^,(Q^(^)) <c,i?-"+("+2)/P vp>o ype[i,^). 

Proof From ((3?22)) and dH]) it foUows that 

11(1 - C)«plL2+4/„(„„+.) < 1(1 - C)^pIIIr"+i < CR-"/'. 
Since C e C^{Qb.{Y)), we then find that for all i? > 0, we have 

(3.27) / _ |7;p|2+4/" < ci?-"-2 vp>0. 
jR"+i\Q„(y) 

Let Ar ^ {X e M"+i : |i;''(X)| > t} and set i? = r-^/". Then 

1^. \ QHiY)\ < / < Cr-^r^ ^ Cr"^ 

■JA^\Qn{Y) 

and n Q/?(y)| < Ci?"+2 Cr-"^. We have proved dS^H). Next, by utilizing 
p.25p instead of (|3.18p and proceeding as in the proof of Lemma 13.31 we obtain 

JQrXY) 

from which (|3.26p follows. The lemma is proved. □ 
Let us fix g e (1, 2±|). By (13:26)1 and ((3T9l) . we have 

By a diagonalization process, we obtain a sequence {p^ij^i tending to zero and a 
function r(-, Y), which we shall call a fundamental matrix for £, such that 

(3.28) TP^{-,Y)^T{-,Y) weakly in Wl'^ {Q R{Y)f Vi? > 0. 
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Then, for any e C;?°(R"+i)^, it follows from jSJ 



= — lim 

fl — ^oo 



(3.29) 



= lim 



Let tJp be the fc-th column of r''(-, y) as before, and let v be the corresponding fc-th 
column of r{-,y). Then, for any g G L^{Qr{Y))'^, ((3:26)) yields 



(3.30) 


/ vg 


= lim 






JQr{Y) 


— >oo 


JQr{Y) 



< Cpi?-"+("+^)/P ||g|| 



iP'(QR(i')) ' 



where p' denotes the conjugate exponent of p e [1. ^^^). Therefore, we obtain 

(3.31) ||r(,y)||^,(Q^(^„ <c,i?-"+("+2)/^' ype[i/-^). 

By a similar reasoning, we also have by (|3.19p 

(3.32) \\DT{;Y)\\^,^Q^^y^^<C,R—'+^-+'yP Vp G [1, ^). 
With the aid of p.27|) . we also obtain 



(3.33) 



|r(-,y) 



li2 + 4/„(R„ + l\Q^(y)) 

Also, by p.22p . Lemma lOl in Appendix, and the assumption ( e C^(Qj^(F)), we 
find (by passing to a subsequence in (|3.28p if necessary) 

(3.34) |||r(,r)|||2„^,^^^(^) <Ci?^-. 



Moreover, arguing as in the proof of Lemma [3.41 and Lemma [3.31 we see that the 
estimates ((XMI) and (|5:M)) imply that for all Y G R"+\ we have 

(3.35) \{x G : |r(A:,r)| >t}\< Ct--^, 

(3.36) I {a: G M"+i : \Dj:{X,Y)\ > t}| < Ct"^. 

Next, we prove in Theorem O Let ^ G Cf{W'^'^) be such that 77 = 1 on 

Qr{Y) for some r > 0. Following the derivation of p.22p . we find 

(3.37) 1(1 - r7)r''(-,r)|||R„+i < C{^) < oo uniformly in p>0. 

Therefore, by passing to a subsequence if necessary, we may assume that there is a 
function r{-,Y) in V2(K"+^)^''^ such that 

(3.38) (1 - r])TP>^{-,Y) f (•,r) "very weakly" in 1/2(R"+^)^''^- 
Note that p.38p in particular implies that 

lim / (i-77)r''-(-,y)</)= / ti;Y)ct) y(f>ec^{w+^f. 

On the other hand, by ((3:281) we find for aU 4> G C;?°(M"+i)^, 
lim / (1 -r;)r''^(-,r)0 = lim / TP" {■ ,Y)(1 - 7])ct) 

= / r(.,r)(i-,?)0= / (i-ry)r(.,r)0. 

JR"+1 JR"+1 
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Therefore, f{-,Y) = (1 - 7?)r(-, y), and thus (l-r;)r(-,r) e y2(K"+i)"^^. Now, 
we win show (l-7?)r(-, Y) e V^-°{W+^)'^'"^. To see this, let v be the fc-th column 
of r(-, Y) and set m = (1 - ri)v. Fix to such that = on (— oo, to] x M". Observe 
that for all T > to, u is a, weak solution in V2{{to, T) x M")^ of the problem 

CU = f + DaOa 

u{to, •) = 0, 

where / := —i]tv + Da'qA"'^ Dpv and := Dp-qA^'^v. Since 77 = 1 on Qr{Y) and 
77 is compactly supported, (|3.34p imphes that /,fif„ G L^{{to,T) x fi)^. Then, by 
Theorem 4.2 in [2D1 §111.4], we have u e V"2^'°((io, x n)'^ . Note that ((337)) implies 
11^11(40, T)xB'» < C'('7)- Since T > io is arbitrary, we have u e V'2^'°((io, 00) x fi)''^. 
Moreover, as in Section [5TT1 if we extend u to the entire R""*"^ by setting it = on 
(-00, to) X R", then we have u e V^'^ 0^'"''^^)'^ ■ We have proved ((2:241) . 

3.4. Continuity of the fundamental matrix. With the aid of (|2.2ip . it follows 
from ((3:29)) and ((334() that T{-,Y) is Holder continuous in R"+i \ {Y}. In fact, 
by the same reasoning, it follows from (|3.3|) and p.20p that on any compact set 
K (£ \ {Y}, the sequence {r''^ (•, F)}^^ is equicontinuous on K. Also, by 

Lemma 12.61 and ((3.27P , we find that there are Ck < 00 and pK > such that 

l|r''(-,y)l|L-(K) < yp<PK for any compact X <E M"+i \ {r} . 

Therefore, we may assume, by passing to a subsequence if necessary, that 

(3.39) r''" (•, Y) ^ r(-, Y) uniformly on K for any compact K d R"+^ \ {Y} . 

We will now show that T{X, •) is also Holder continuous in \ {X}. Denote 
by *r'^(-,X) the averaged fundamental matrix associated to *£. More precisely, for 
1 < fc < iV and X = {t, x) G R"+^ fixed, let w^, — Wa--x,k be the unique weak 
solution in y2^'°((— cx3, <o) x R")^ of the backward problem 




where to > t + tr^ is fixed but arbitrary. Then, as before, we may extend Wa to 
the entire R"+^ by setting Wa = on (t + a'^,00) x R" so that belongs to 
V2 ' (R"^^) and satisfies for all Si < t the identity 



'2 

(3.40) / wl<t>\sr,-)+ Tj wl<t^\+ ri 'AfDpwlD^^^ = i <\>^ 
for all G C~(M"+i)^. We define the averaged fundamental matrix *r''(-,X) = 

(*rjfe(-,x))^Vi for*/: by 

*rj,(,x)^7«i = <^_,. 

By a similar argument as appears in the previous subsection, we obtain a sequence 
{dvS^^x tending to and a function 'r(-,X), which we shall call a fundamental 
matrix for *£, such that for some g > 1, 

(3.41) T'"''(-,X) ^ *r(-,X) weakly in VKi''?(Qfl(X))^^^ Vi? > 0. 
Moreover, as in ((3.29p . *r(-,X) satisfies 

T,fc(-,x)0j + / X%*r,fe(.,x)D„<^' = 0'=(X) 



18 



S. CHO, H. DONG, AND S. KIM 



for all 4> G C^(M"+^)^. Also, by a similar reasoning, we find that *r(-,X) satisfies 
all the properties corresponding to (I3.3ip - (|3.36|) and is Holder continuous in R"+^ \ 
{X}. Therefore, the following lemma will prove the claim that r(X, •) is Holder 
continuous in \ {X} as well as the estimates (IZ;28)) - ((237)) . 

Lemma 3.5. LetT{-,Y) and^{-,X) he fundamental matrices of L and^C, respec- 
tively, as constructed above. Then 

(3.42) T{X, Y) = T{Y, X f VX Y. 

Proof Fix X,Y e R"+i with X ^ Y. From the constructions of r''(-,r) and 
*r'^(-, X), it is not hard to verify that 

(3.43) / *r^,(.,x) = / rf,(,y). 

JQ7(Y) JQt(X) 
Let {Pfj,}^^i and {a-^}^^^^ be approximating sequences for T{-,Y) and *r{-,X), 
respectively, and denote 



From the continuity of r^^(-, Y), it follows that 



Ji^^M. = £^ / , C (•, r) = rf,r [x, y), 



'QZAX) 

and thus, by ([3?39| . we obtain 



lim lim lim V^- {X ,Y) = T ik{X ,Y) . 

On the other hand, (I3.4ip yields that 

lim g';i - lim / Kii'^X) = / ^Tkii^X), 

and thus, it follows from the continuity of ^Tkii'TX) that 

lim lim g';l^ lim / H:u{-,x) ^H:ki{Y,x). 

/J. — *-oo f — >oo ^ /J. — >oo J Q- (^Y) 

We have thus shown that Tik{X, Y) = Tfc; (F, X) for aU X ^Y. □ 

So far, we have seen that there is a sequence {Pi_i}^^i tending to such that 
rP''(., Y) r(-, Y) in \ {y}. However, by ^Ml, we obtain for X ^ F, 

rf,(x,r) ^ hni^/ rf,(-,n = / ^rO.^) 



Tkli;X)^4 TikiX,-). 

iQp{y) ■JQpiY) 

Therefore, we have the following representation of the averaged fundamental matrix: 

(3.44) r''(x, Y) = -l r(x, z) dz yx ^ y. 

jQpiY) 

In particular, by using the continuity of T{X, •), we obtain 

(3.45) iimr''(x,y) = r(x,r) yx^Y, 
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and thus, from (I3.17P that 

(3.46) \T{X,Y)\ <C\X -Y\-'' VX^Y. 
Furthermore, due to p.2p . it also follows that 

(3.47) r(-,-,s,?/) = in {-oo, s)xW\ 

3.5. Representation formulas. We shall prove the identity (|2.22p . For a given 
/ G C;?°(R"+i)^, fix to such that / = on (-oo,to] x R". Let u be the unique 
weak solution in V2'^'°((io, oo) x R")^ of the problem 

Cu^f 

By the uniqueness, u does not depend on the choice of to and we may extend u to 
K"+i by setting u = on (-oo,io) x E" so that u e V^'°{W+^)^. Let 'T''{-,X) 
be the averaged fundamental matrix for Then, as in (13. 8|) . we have 

(3.48) / *r,^(,x)r = / 

As in ()3.12p . the property (PH) implies that u is continuous. Therefore, by taking 
the limit ^ oo in p.48p and then using (I3.42p . we find that 

rM{x,-)f ^u'^ix) 

or equivalently, in terms of matrix multiplication 

u{X)= [ T{X,Y)f{Y)dY. 

We have thus proved ([2?22ll . 

Next, we shall prove the identity fTI^ . Let g e ^^(K")^ be given and let 
u be the weak solution in V"2^'''((s, oo) x R")^ of the problem ((2:26)) . Fix X = 
(i, x) G R"+^ with t > s and let Wa = Wa;X,k be the fc-th column of the averaged 
fundamental matrix of the transpose operator *£ with a pole at X. By using ()3.40p 
and (|2.8p together with standard approximation techniques (see e.g., [20l §in.2]), 
we find that for Ui, sufficiently small, 

(3.49) / u""^ ( w,,is,y)-9{y)dy. 

Let us assume for the moment that g is compactly supported. Then, by estimates 
similar to and (piiS)) . 

lim / w„^{s,y) ■ g{y)dy ^ I w{s,y) ■ g{y) dy, 

where w is the fc-th column of *r{-,X). Also, by the property (PH), we find that 
u is continuous at X. Therefore, by taking the limit — > oo in (|3.49p . we obtain 

(3.50) u\X)= [ *r,k{s,y,t,x)g\y)dy^ [ rk^{t,x,s,y)g'iy)dy. 



We have thus derived ()2.25p under the additional assumption that g is compactly 
supported. For general g e i^(R")^, let {gm}m=i be a sequence in L^(R")^ 
such that limm_^oo lISm — g|li2(K„) = 0. Let 

Urn be the unique weak solution 
in V2^'°((s,oo) X R")^ of the problem ((2:26ll with g replaced by g^. Then, by 
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the energy inequality, we find that hnim^oo Wum — M|||(s,i)xR" = 0. Therefore, by 
Lemma [2.61 we conclude limm^oo \urn{X) — u{X)\ = 0. On the other hand, by 
([2:24| and ((XSi)) . we have ||r(i, x, s, •)|li2(R„) < C(i-s)-"/2, and thus 

lim / T{t,x,s,y)g^{y)dy ^ I T{t,x, s,y)g{y) dy. 

This completes the proof of the representation formula ()2.25p 

Next, forg e L^{W)'^ 
of the backward problem 

(3.51) 



Next, for g e L^(R")^, let u be the unique weak solution in F2^'°((-oo, t) xR")^ 



*£n = 
u{t,-) 

Then, by a similar reasoning as above, u has the following representation 
(3.52) u{s,y)^ [ *r{s,y,t,x)g{x)dx. 



3.6. Uniqueness of fundamental matrix. We have already shown that T{X, Y) 
is continuous in {{X, Y) e M"+i x K"+^ -.X^Y} and satisfies T{t, x, s,y) = for 
t < s. Also, we have seen that T{X, •) is locally integrable in R"+i for aU X e R"+i 
and that for all / S C;?°(K"+^)^, the function u given by iT22\i belongs to the 
space y2^'''(R"+^)^ and satisfies Cu — f in the sense of (j2.5p . Suppose F{X,Y) is 
another fundamental matrix satisfying all the properties stated above. Then, for a 
given / G C;?°(R"+i)^, let u be defined as in (|2:^ and set 



u{X) j^^^f{X,Y)f{Y)dY. 



Since T{t^x,s,y) = T{t,x,s,y) = for t < s and / is compactly supported in 
we find it = m = on (— oo, to] x M" for some to. Therefore, v := u — u is in 
the space V2'°ii'to, oo) x K")^ and satisfies £i; = in (to, oo) x R" and v{tQ, •) = 0. 
By the uniqueness (see [101 §111.3]), we must have v = 0. Therefore, we find 



/ {r~f){X,Y)f{Y)dY^Q \ffeC^{l 



and thus it follows that r(X, •) = f{X,-) in R"+i \ {X}. Since X e W'+^ 
arbitrary, we have proved that T = r in {{X,Y) e R"+i x R"+i : X ^Y}. 



3.7. Conclusion. We completed the proof of Theorem l2.7l in the case when = R'^ 
and Rc ^ oo, modulo the identity ([2?271) . We defer the proof of ([2?27| to SectionO 

4. Proof of Theorem 12.71 General cases 



In this section, we prove Theorem 12.71 under a general assumption that is 
an arbitrary open connected set in R" and Rc £ (0,oo]. Since we allow the case 
n ^ R", we shall write G{X,Y) instead of T{X,Y) for the Green's matrix in 
U ^ Wx CI to avoid confusion. To construct Green's matrix in U, we need to adjust 
arguments from the previous section a little bit. 

Throughout this section, we shall use the notations Up{X) ~ Q^{X) OU, 
U+{X) = Q+{X) nU, and Up{X) = Qp{X) nU. Also, recall that we have de- 
fined dx = dist{X,dU) = dist(a;,9ri) and dx — niin{dx,Rc)- As in the previous 
section, we employ the letter C to denote a constant depending on n, N, A, A, /io, 
Co, and sometimes on an exponent p characterizing Lebesgue classes. 
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4.1. Averaged Green's matrix. Let Y ^ {s,y) E U and 1 < fc < be fixed. 

For each p > 0, fix sq G (— oo, s — p^). We consider the problem 



(4.1) 



m(so,-) = 0. 



By Lemma |2. 11 there is a unique weak solution Vp ~ Vp-Y,k in V2'^'*'((so, oo) x 
of the problem ()4.ip . As in Section [321 we may extend Vp to entire U by setting 

(4.2) = on (-00, s - p"^) X il 

so that Vp G V2'^{U)^ and satisfies for all ti > s the identity 

(4.3) / vl(p\t,,-)- r / r / AfDpvlD^^^ = i 0'= 

Jf2 J-ooJn J-QoJQ. JUp(Y) 

for all € C^p{U)^ . We then define the averaged Green's matrix G''{-,Y) = 

Next, for each / G L'^{U)^, fix such that / = on [to,oo) x fl and let u 
be the unique weak solution in V2'^'"((— oo, io) x il)^ of the problem (|3.4p . Again, 
we may extend it to Z/^ by setting m = on (^o, oo) x il so that u G V2'^{U)^ and 
satisfies for all ii the identity 



(4.4) /^,Y(ii,-)+/ 'AfDpuW^cj,^ ^ fcj,^ 
Jn Jti Jn Jti Jn Jti Jn 

for aU G C^p(Zi). 

We note that Lemma [01 remains valid with il, U, and Z^^" in place of M", K"+^, 
and Qp , respectively. By following a similar argument as in Section 13. 2i we find 
that if / is supported in Q'^{Xq), where Xq G U and R < dx„, then we have 

II«IIl=o(q+^^(x„)) < CR'-^''^'^^^fh.(QUXoy) yp>in + 2)/2. 
Consequently, as in (|3.16p . if we have p < R/A and Qp{Y) C Q^^^{Xo), then 

||^plL,(Q+(x„)) < Vg G [1, (n + 2)/n). 

Therefore, by following the proof of Lemma 13.21 we conclude that 

(4.5) \GP{X,Y)\<C\X -Ylp'^ yp<\X~Y\p/3 if |X - y|p < (iy/2. 

4.2. Construction of the Green's matrix. As in p.20p . the estimate (|4.5p vields 

(4.6) \lGP{-,Y)l\l^^^^^^<CR-" Vp>0 Vi?<Jy/2, 

which in turn gives the following weak-L("+^^/'^"+^^ estimate as in p.lSp : 

(4.7) \{X eU ■.\D^G''iX,Y)\> t}\ <Ct-^ Vt > (Jy/2)-". 
Then, as in (|3.19[) . we have for any p > 0, 

(4.8) pG''(.,r)||^,(Q^(^))<Cpi?— VRKdy VpG[l,^). 
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Also, as in (|3.27|) we have 



(4.9) /_ |G''(-,r)|'+^/" < Ci?-"-2 y^>o Vi?<dy/2, 



JU\Q^{Y) 



which gives the following weak-L'^"+^^/" estimate as in p.25p : 

(4.10) \{X^U: \Gf{X,Y)\ > t}\ < C7r-("+')/" Vr > (Jy/2) 



— n 



Then, as in ()3.26p . we have for any p > 

(4.11) ||G''(,y)||^,(Q^(y))<Cpi?-"+("+2)/P vi?<dV Vpe [l,(n + 2)/n). 
Fix q e (1, ^) and observe that (|48l) and (|4?TT1) imply 

(4.12) ||G''(-r)||;^^i,,(Q^. ^Y)) ^ C{dY) < oo uniformly in p > 0. 
Also, fix a cut-off function 

(4.13) Cy eC^{Qs,/^{Y)) satisfying (y = I on Q^^ /,{Y). 
Then, as in (|3.37p . we find 

(4.14) 1(1 - CY)Gf{;Y)\lu < C{Cy) < oo uniformly in p > 0. 

Therefore, we conclude from (|4.12p and (|4.14p that there exist a sequence {Pi_i}^^i 

tending to and functions G{-,Y) e W^'i{Qa^{Y))^''^ and G{-,Y) G V2{U)^''^ 
such that as /X — > oo, 

(4.15) G''-(-,r)-G(.,r) weakly in W^i^«(g,-^(r))^><^ 

(4.16) {1~(:y)GP'^{-,Y) ^G{-,Y) "very weakly" in F2(Z^)^'' ^■ 

Since we must have G{-,Y) = G{-,Y) on Qjy(F) \ Qjy/2(i^), we may extend 
G(-,y) to the entire by setting G{-,Y) = G{-,Y) onU\Q^^{Y). We shall caU 
this extended function a Green's matrix for C'vaU with a pole at Y and still denote 
itbyG(-,y). 

Now, we prove the identity (|2.23p in Theorem 12.71 Write 4> = ^0 + (1 ^ '7)0> 
where r/ e C^{Qa^{Y)) satisfying 77 = 1 on QdV/2(>"). Then, by gl3|), gT5l) . and 
(|4.16p , we obtain 

(t)^{Y)= lim / 77(/)'=+ lim / (1 - ?7)0'' 




+ lim 



/ (-G'^i:{-,Y)D,{{l-riW)+Al^D,G% 

J lA 



(•,r)i^„((i-7?)^')) 
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By proceeding similarly as in Section l3.3|, we derive the following estimates which 
correspond to ()3.31|) - ()3.36p : 

(4.17) ||G(.,y)||^,(Q^(y)) <C,i?-"+("+2)/'' VRKdy ype[l,^), 

(4.18) pG(.,r)||^,(Q^(^)) <C,i?— yRKdy Vpe[l,^), 

(4.19) l|G(-, r)||i.+4/„(;,\Q^(y)) < Ci?-" Vi? < Jy/2, 

(4.20) |||G(.,r)|||^^^^(^j <Ci?-" Vi?<dV/2, 

(4.21) |{X e : |G(X,r)| >r}| <Cr~'^ Vt > (dV/2)-", 

(4.22) \{X eU ■.\D^G{X,Y)\>t}\<Ct-'^ Vt > (Jy/2)-". 

Next, we prove ((2:24)) in_Theorem [2Jl Let ry € (W) such that 77 = 1 on (5^(5^) 
for some r < dy- Since dy < dy, we may assume that r < dy. Then, as in (j3.38p . 
we may assume that there is a function G{-,Y) in V2{U)^^^ such that 

(4.23) (1 -77)G''^(-,y) ^ G(-,y) "very weakly" in \/2(Z^)^''^- 
Note that (|4.23p in particular implies that 

hm / {l-Tj)GP-{; ¥)<(>= f G{;Y)4> Vcj^eCnU). 

On the other hand, by (jil^ and dHH), we find that for aU cj) e C^{U), 

lim / (1 - ry)G''^(-,y)0 

= lim ( / G''-(-,y)(l-77)(/.+ / (l-Cy)G''-(-,y)(l-?7)0) 
= / G{-,Y){1 - 77)0 + / G(-, y)(l - ry)0 - / (1 - ry)G(-, y)0. 

"'QdV(^) ''U\Qi^(Y) Ju 

Therefore, we must have G(-,y) = (1 - r/)G(-, y), and thus (l-r/)G(-,y) G V2(Z^). 
Then, by following a similar argument as in Section [3731 we find that (1 — ri)G{-, Y) 
belongs to V^'^W^ which proves ([2^. 



4.3. Continuity of the Green's matrix. Proceeding as in Section [331 we find 
that G(-, y) is Holder continuous in U \ {Y}. Also, we may assume that 

(4.24) Gf^ (•, y) ^ G(-, y) uniformly on K, for any compact K (^U\{Y}. 

Denote by HS''{-,X) the averaged Green's matrix associated to *£ with a pole at 
X £U. More precisely, for 1 < fc < iV and X — {t,x) & U fixed, let i^o- = Wa-x,k 
be the unique weak solution in V2 ' ((~oo,to) x il)^ of the backward problem 

I " |w+(x)| 
1 M(io,-) = 0, 

where to > ^ + "'^ is fixed but arbitrary. Then, as before, we may extend Wg^ to the 
entire U by setting = on (< + cr^, 00) x fl so that lOo- belongs to y2^'"(Z^) and 
satisfies for all si < t the identity 

(4.25) / <</.''(si,.)+ r/ <0j 
Jo Jsi Jn Jsi Jn " Jui(x) 



r[ 'AfDpwiD,^^^ ^ / 
Jsi Jn Jui 
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for all e C^^{U)^. The averaged Green's matrix *G'"(-,X) = {^G''^ki-^^))lk=i 
for *£ in U is then defined by 

By a similar argument as appears in the previous subsection, we obtain a sequence 
{cr^i;^! tending to and functions G Wl^'i{Q^^{X))^'^^ for some g > 1 

and ^G{-,X) e ^iU)"""^ such that as -> oo 

(4.26) *G^-(.,X) - *G(.,X) weakly in W^^'^iQ^^iX))^-^ , 

(4.27) (1 - Cx)*G^-(-,^) ^ *G(-,^) "very weakly" in V^iUf'', 

where is given similarly as in (|4.13|) . Since we must have 'G(-,X) = 'G(-,X) on 
Qdx i-^) \ Qdx/2{^)j extend *G{-,X) to the entire U by setting *G{-,X) = 

G{-, X) onU\Q^^{X). We shall call this extended function a Green's matrix for *£ 
in U with a pole at X and still denote it by *G(-, X). Then, by the same reasoning 
as in the previous subsection, we find that *G(-, X) is Holder continuous mU\ {X} 
and it satisfies 



(^'G,k{-,X)cl)l + 'A^fDt,'G,k{-,X)D^q^') = /(X) V0 e C^U) 



N 



and all the properties corresponding to (|4.17p - (|4.22p . Therefore, the estimates 
(|2.28p ~ (|2.37p will follow once we establish the following identity, the proof of which 
is essentially given in Lemma 13.51 

(4.28) G{X, Y) = *G(y, Xf yX,Y eU, X^ Y. 

Moreover, (14.281) also implies that G(X, •) is Holder continuous in U \ {X}. Fur- 
thermore, as in (I3.44P and ()3.45p we have 

G''{XX) = l G{X,Z)dZ \/X,YeU, X^Y and 

Ju-{Y) 

\imGP{X,Y) = G{X,Y) yX,Y eU, X ^ Y. 
Then, from and ((i?^ , we find 

(4.29) |G(X,y)| < C|X -r|p" if < |X - r|p < imax(Jx,rfy) 
and also by ()4.2p . we obtain 

(4.30) G(-,-,s,y) = in {~oo,s)xn. 

4.4. Proof of the identity ([2^. Let g G L'^{n)^ be given and let u be the 

unique weak solution in y2^'°((s,oo) x Vl)^ of the problem (j2.26p . By proceeding 
similarly as in Section [3.51 we derive the following representation for u{t,x): 



(4.31) u{t,x)= G{t,x,s,y)g{y)dy. 

Jn 

Similarly, for g e L^(f2)''^, let u be the unique weak solution in V2'^''^((— oo, t) x 17)^ 
of the backward problem p.5ip . As in (|3.52p . u has the following representation: 

(4.32) uis,y)= [ 'G{s,y,t,x)g{x)dx. 

We need the following lemmas to prove (|2.27p of Theorem 12.71 
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Lemma 4.1. Assume that g G L'^{fl)^ is supported in a closed set F C ^l. Let u 
be the weak solution in V2'^i{Si'^) ^ i^)^ of the problem p.26p . Then, 

(4.33) /" |tt(t,-)|^ < e-^'i'''*(-^'-^)'/(*-^) /" ^Ecn, 

Je J f 

where dist(i?, F) = inf {|a; - y\ : x e E,y e F} and c = \/{2IS?). 

Proof. We may assume that dist(i?, F) > 0; otherwise (|4.33p is an immediate con- 
sequence of the energy inequahty. Let V' be a bounded Lipschitz function on f2 
satisfying \Dip\ < 7 a.e. for some 7 > to be fixed later. Denote 

I{t) / e^'^\u{t,-)\^ yt>s. 



By following '20', §111.4], it is not hard to see that / is absolutely continuous on 
[s,oo) and that I'{t) satisfies for a.e. t>s 



I'{t) = -2 J^^ A'tj'Dpu^(t, ■)D^{e^^uHt, •)) 

= -2 [ e^'^A%PD0U^t,-)D^u\t,-)~A [ e^^ A^f D^u^t, ■)D„i; u\t, ■) 
Jn Jn 

< -2A / e^'l'\Du{t,-)\^ +Ak-i [ e"^ \Du{t, ■)] e"^ \u{t, ■)] 
Jn Jn 

< (2AVA)72 / e^'^ \u{t,-)\' = {2K^ / \)^^I{t). 
Jn 

The above differential inequality yields 

(4.34) I{t) < e(2A'/^)^'(*"^)||e'^g||i.(j,). 

Note that since F is closed, the function dist(a;,F) = min{|a; ^ y\ ■ y E F} is a 
Lipschitz function with Lipschitz constant 1 and dist(£',F) — inf^jg^ dist(a;, F). 
Therefore, if we set ipix) — 7 min{dist(x, F), dist(£^, F)}, then by (|4.34p . we find 



/ |M(t,-)|' <exp{(2AVA)72(t-s)-27dist(F,F)} / \g\\ 
Je J f 

The lemma follows if we set 7 = dist(£;, F)/{(2AVA)(t - s)}. □ 



Lemma 4.2. Assume that the operator C satisfies the property (PH). Let u be the 
weak solution in V2 '^{{s, 00) x Vi)^ of the problem (j2.26p with g £ L'^{ fl)^ . Then, 



|M(i,a;)| < C||g||i„o(o) tfVt~s< dx, 
where C = C{n, N, X, A, fio,Co) > 0. 



Proof. Let p := \/t — s and assume p < dx. Denote ^0 = Bp{x) and Ak = 
{y £ n : 2*^^^ p <\y — x\ < 2*''p} for k > 1. Since g is compactly supported in il, 

we see that g = X]fc=o 91-41- f*-"^ some fco < 00. Denote 



Uk{t,x)= G{t,x,s,y)g{y)dy. 
J A. 



Then, it follows from (|4.3ip that u — X^fclo and that each Uk is a unique weak 
solution in V^^'°((s,oo) x il)^ of the problem (|2.26p with glA^ hi place of g. If we 
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apply Lemma |4?T] to Uk with E = Bp{x) and F — Ak, we obtain 

\uk{r.y)\^ dy < Ce-^^^^^p^ \\9\\l^<n) Vr e (s,t). 



Therefore, by Lemma 12.61 we estimate 

\u(t,x)\ < ^ \uk{t,x)\ < C^e"=2'=2'="/2 < C^||9lli=o(n) ■ 

k=0 fc=l 

The lemma is proved. □ 

Lemma 4.3. Let xq E il and R < ^ min(dist(a::o, (?f2), i?c) be given. Assume that 
T] £ C'^{B2r{xq)) satisfies < r/ < 1, ry = 1 m Bji{xn), and \Dr]\ < C/R. Then, 

lim / G{t,x,s,y)rj{y)dy ^ I ^x e Br{xo), 

where G(t,x, s,y) is the Green's matrix of L as constructed above and I is the N 
by N identity matrix. 

Proof. First, we claim that the following identity holds for all si < t: 

(4.35) </)'=(X) = / *G,fc(si,-,i,x)0*(si,-)+ / f 'G,k{-,X)cj)l 

+ f f *A^fDp'G,k{;X)D^^' V0 G G^piU)"" . 
JsiJn 

Indeed, by (|4.25p . we have for all si < i, 

(4.36) / <jy^ ^ hG2{si,-.t,xW{s^,-)+ rhG2{;X)<p\ 

'A^fDp'G^^i: {■,X)D^f V0 e (Zi)^. 

Note that by estimates similar to (|4.24l) and (|4.29[) . we have 

lim / *G'^^{si,-,t,x)(t)''{si,-) ^ / ^G^k{sl,■,t,x)(j)\sl,■). 

Then, (|4.35p follows from a similar argument as used in the proof of (|2.23p in 
Section I4?2l and ((430)) . If we set (f){Y) = rj(y)ei and si ^ s < t in ((435| . then 

4i= f *Gik{s,;t,x)r,+ f / *A"%*G,fe(-,X)7?,r; 
Jn J s Jn 

(4.37) = / Gki{t,x,s,-)r]+ f f 'A"" Dp'G,k{-, X)D^r, =: I + 11. 

Jn Js Jn 

By using Holder's inequality and (|4.20p . we estimate 

|//| < Ci?"/2-i(i-s)i/2 ( f[ \D'G{-,Xo)i\ 



\J sJ BR{xo)\Bfi/2{xo) J 

<CR-\t~s)^/^ ifVt^<R. 
Therefore, the lemma follows by taking the limit i to s in (|4.37p . □ 
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Now, we shall prove ([2^ of Theorem [2Jl Assume that g E L'^{VL)^ and 
g is continuous at x^j. Let u be the weak solution in V'2^''^((s, cx)) x 51)^ of the 
problem (I2.26p . For a given e > 0, choose R < ^ min(dist(xo, c)ri), i?c) such that 
\g{x) — g{xo)\ < e/2C for all x satisfying |a; — xo| < 2R, where C is the constant 
as appears in Lemma 14.21 Let rj be given as in Lemma 14.31 and let Mq, u^, and 



Moo be the weak solutions in V^^'''((s, oo) x fi)^ of the problem (|2.26p with 7]g{xo), 
vid ~ g{^o))j a-nd (1 — rj)g in place of g, respectively. Then, by the uniqueness, we 
have M = Mo + Me + Uoo- First, note that by (|4.3ip . Uq has the representation 

(4.38) Uo{t, x) = (^J G{t, X, s, y)ri{y) dy^ g{xo). 

Next, we apply Lemma |4. II to Mqo with E = B]i/2{xo) a-^id _F = f2 \ Be.{xq) to find 
|«co(r,y)|' dy < e-^^-Vf-^^) Vr > s. 

Therefore, if p := \/t — s/2 < i?/4, then by Lemma [2?6l we find 

(4.39) |Moo(i,x)| <Cp-"/2e-^(«/'')'||g||^.(f,) Vx e Bp{xo). 
Finally, we estimate by using Lemma [ 



(4.40) \u^{t,x)\<e/2 iiVt^Kdx- 

Combining (|4.38p . (|4.39p . and (|4.40p . we see that if \/t — s is sufficiently small, then 
for all x e Bp{xo), we have \u{t, x) — g{xQ)\ < e. This completes the proof. 

4.5. Conclusion. The proof of representation formulas (|2.22p and (|2.25p given in 
Section 13.51 as well as the proof of the uniqueness given in Section 13.61 also works 
for general domains O and Rc € (0,oo]. Therefore, Theorem 12.71 is now proved. 

5. Proof of Theorem 12.111 

5.1. Proof of the Gaussian bound ()2.40p . Here, we consider the case Rc — oo 
and prove (I2.40p . In [l4j, by following methods of Davies ;4^ and Fabes-Stroock 
Hofmann and Kim derived the upper Gaussian bound of Aronson 1 under 
a further qualitative assumption that the coefficients of C are smooth. For the 
reader's convenience, we reproduce their argument here, dropping the technical 
assumption that the coefficients are smooth. 

Let -0 be a bounded Lipschitz function on M" satisfying \D4>\ < 7 a.e. for some 
7 > to be chosen later. For t > s, we define an operator Pf_^t on _L'^(M")^ as 
follows. For a given / e L'^{W')^, let u be the weak solution in F2^'°((s, oo) x M")^ 
of the problem 

^rc 1^ ( Cu = 

Then, we define Pf_^tf{x) := e'l'^''^u{t,x). Note that it follows from ^(n^ that 
(5.2) pLJ{x) = e'^(^) / r(t, X, s, y)e~^^y^f{y) dy V/ G i^^^")^. 

We denote I{t) := ||e'^M(t, •)||i2(R„) for aU t > s. Then, as in the proof of 
Lemma |4T| we find that I is absolutely continuous and satisfies for a.e. t > s, 

I'{t) < (2AVA)7'/(i). 
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The above differential inequality with the initial condition /(s) = II/IIl2(r,i) yields 

< ,(2AVAh^(*-) yt>s. 
Since I{t) = \\Pf_^tf\W2(^f^r,-^ for t > s, we have derived 

(5.3) ll^'.tt/llL^(E") <e'"''(*-^Ml/llL2(M.) yt>s, 

where v = A^/A. By Lemma we estimate 



C 



Hence, by using (|5.3p we find 



We have thus derived the following L°° estimate for Pf_^). : 

(5.4) ||P,t,/|U»(E") <C(t-.)-"/^e^v^+''^^(*-^'||/llL.(M.) Vt>s. 

We also define the operator Qf_^s '^^ i^(M")^ for s < i by setting Qf^^g(y) = 
e~'^(y^v{s, y), where v is the weak solution in ¥2 '^ {{—oo, t) x M")^ of the backward 
problem 

Then, by we find that 

Qf_.g(y) = / *r(s, y, t, a;)e'^(-)g(:r) dx Vg E L'{WY- 

By a similar calculation that leads to (|5.4p . we derive 

(5.6) ||Qt,g|U^(K.)<C(t~s)-"/4e'^^^+''^'(*-^)||g||^.(„,.) < t. 
Notice that it follows from and that (cf. ([SI]) and ([Qg]) ) 

/ (/'.tt/)-9= / /•(Qt.s). 

Therefore, by duality, (|5.6p implies that for any / e i^(M")^, we have 

(5.7) WPtJWmu'^) < c{t - ,)""/4g7y*^W(*-) ||/||^,(^„^ vt > s. 
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Now, set r — (s + i)/2 and observe that by the uniqueness, we have 

PLtf ^ PLt{Pt-^rf) V/eL-(M«)^. 

Then, by noting that t — r = r — s — {t — s) /2, we obtain from (|5.4p and (|5.7p that 
for any / e if (M")^, we have 

For fixed x,y ^ M" with x ^ y, the above estimate and (|5.2p imply, by duahty, 

(5.8) e^(-)-^fe) \Y{t, X, s, < C{t - s)-"/^ eV2(t-^) > s. 

Define ^/ig on [0, oo) by setting 

f V'o('') = J' ''iir <\x-y\ 

\ Mr) ^\x-y\ if r > jx-yj. 

Let il>{z) :— 7V'o(|2 — y\) where 7 = \x — y\ /2h'{t — s). Then, is a bounded 
Lipschitz function satisfying \Dip\ < 7 a.e., and thus (|5.8p yields 

|r(t, x, s, y)|„p < C{t - exp{e/A/2i. - ^VM, 

where C = 1^^ ~ ?/l l\Jt — s. Let Co = Co(t^) = Co(A^/A) be chosen so that 

exp(,e/\/2i^ - eV4i^) < C'o exp(-eV8i^) G [0, 00). 
If we set K — Xj'Ssv = A/8A^, then we obtain 

|r(i, X, s, < C{t - exp {-kIx - yfl{t - s)} , 

where C = C(n, TV, A, A, /xq, C'o) > 0. We have proved ((2^ . 

5.2. Proof of (p:iT1) and 1^:^ . We begin by proving the identity (P^ITjl . First 
note that (H^, ^T^ . and (P?^ imply that both |r(t, r, •)| and |r(r, -,5,^)1 
belong to for all x,y ^ R" and r € (s,i). Then, by the uniqueness (cf. 

and Fubini's theorem, we have for aU g e i^(R")^, 

/ ^{t,x,s,y)g{y)dy = T{t,x,r,z)( r{r, z, s,y)g{y) dy] dz 

Since t/ G L^(R")^ is arbitrary, we conclude ^A^. 

We now turn to the proof of the pointwise bound (|2.39p . By following the proof 
of (|2.40p in Section 15.11 it is routine to check 

(5.9) |r(t,a;,s,y)|„p < C(t-s)-"/2e-''l"-''l'/(*-'^) ifO <t-s<Rl. 
Next, recall the identity 

(5.10) <^{t + s,x-y)= i <^>{t,x- z)^{s,z-y)dz Vx,?/eM" Vs,t>0, 

where $(t, x) = (47rt)""/2 exp{-|x| V4t}. Note that implies that 

(5.11) \T{t,x,s,y)\<C{n/Kr'^^{{t-s)/^K,x-y) iiO<t^s<Rl. 
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Let I be the largest integer that is strictly less than [t — s)/E^. Denote tj = s + jR^ 
for j — 0,...,£ and ti+i = t so that to = s and ti < t ^ ti+i. Note that 
|ij < Rl for all j 1, . . . ,^+1. Then, by (gH]), (j^T^ . and we have 

\T{t,zt+i,s,zo)\ < ■■■ T[\r{t-i,Zj,tj^i,Zj^i)\ dzi--- dzf: 

< {C(7r/^)"/2}^+i$((i „ s)/4aj, z,+i - zo) 

Therefore, l|2.39p follows if we set zq = y, z^+i = x, and 7 = ln{C(7r/K)"/^}. 

5.3. Proof of the identity ^L^. First, note that as in (HSU), the following 
identity holds for all si < i: 

(5.12) 4>\x)^ I 'T,k{su-:t,xW{su-) + fjjr,k{-,xw^ 

Now, let C, e C^(R") be a cut-off function satisfying 

C=linBi(0), C = outside ^2(0), < C < 1, \DC\<2. 
If we set = C{{y - x)/R)ei and si = s < t in (|5.12p . then 

(5.13) 4i= / *rzfc(s,-,i,x)C((--x)/i?)+ // Uf/i?^*r,fc(-,x)i^,0' 

= 1 + 11. 

By (|3.42p . (|2.39|) . and the dominated convergence theorem, we have 



lim /= lim / Tkiit, x, s,y)C{{y ~ x)/ R) dy = / Tki{t, x, s,y) dy. 
On the other hand, using Holder's inequality we estimate // by 

\ii\<CR''^^~\t-s)^^^ ( ff \D*r{-,x)\'''] . 

\JsJB2R{x)\BRix) J 

Next, let 7] £ C'^{B^fi{x)) be such that 77 = 1 in B2b.{x)\Bb.{x), 77 = in B]i/2{x), 
and \Drj\ < C/R. From ([QO]) and (g^ll), it follows that limstt r/*r(s, •, t, a;) = 0. 
Note that T(-,X) satisfies *Cu = in {{s,y) G ]R"+^ : s < i}. Therefore, as it is 
done in the proof of p.2ip . we derive 

\D*r{-,X)\^ <CR-^ f f \*r{-,x)\\ 

I sJ B2r{x)\Br(x) J sJ B3r(x)\Br,2(x) 

If i? > then by ([428)1 and ((239)) . we estimate (cf. ([3^) 

Therefore, we find that // is bounded by 

\II\ < CR-\t - s)l/2g7(t-«)/i?^ if ^ > ^t^^ 

and thus, we obtain ()2.42p by taking i? 00 in (|5.13p . 



|*r(-,x)|' < Ce2^(*-^)/«-i?2- 
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6. Appendix 

Lemma 6.1. Let {wfc}^]^ be a sequence in V2{li) (resp. V2{1^) ), where Li — M.x. Q, 
and f2 is an open connected set in M". If sup f. |||ufc|||w — M < oo, then there exist 
a subsequence {uk^}'j^i C {uk}'^i and u G ¥2(14) (resp. V2{U)) with Wufu < M 
such that Ukj ^ u ''very weakly" in V2(U). 

Proof. Let Ut '.^ U f] {~T,T) x K", where T > 0. We claim that for all T > 0, 
there exist a subsequence {ufe^}°^i C {ufe}^^ and u £ V2{Ut) with |||it|||i/T, < M 
such that u weakly in W2''^{Ut). It will be clear from the proof below that if 

{ufc}^-^ is a sequence in ^2(^)1 then we also have u e V2(^t)- Therefore, once we 
prove the claim, the lemma will follow from a standard diagonalization argument. 

Now, let us prove the claim. Notice that ||ufe|lvi/i'"(z^y) ^ C'lll^felllwT — ^^'^ fo^' 
some C — C{n,N,T) < 00. Therefore, there exist a subsequence {ukj}j^i C 
{wfelfe^i and u e W2'^{Ut) such that such that Uk^ u weakly in WI'^^IAt). For 
simplicity, let us relabel the subsequence and assume that Uj ^ u in W2 ''^{Ut). We 
only need to show that u £ V2{Ut). Let t £ {—T,T) be fixed and choose ft, > so 
small that {t — h,t + h) C {—T, T). Since Uj u weakly in L'^{Ut), we also have 
Uj u weakly in — h,t + h) x fl), and thus 

rt+h p pt+h ^ 

\uj{t,x)\'' dxdr < M'\ 

n 

Therefore, by taking h ^ 0, we obtain ess suptg(_2^ -jUla < Af^, which 

together with u £ 1^2^'° (ZYt) implies that u £ V2{Ut)- It is clear that < 
supj, |||itfe|||;^7, < M. We have proved the claim, and thus the lemma. □ 

Lemma 6.2. Assume the coefficients A"^ of the operator C in (|2.2p are indepen- 
dent of X (i.e., A"^ — A°'^{t)) and satisfy the conditions (12. 3|) and ()2.4p . Then, 
there exists a constant C = C{n, N, X, A) > such that all weak solutions u of 
Cu — in Q^{Xq) satisfy 

|2 ^ ^ /p\"+2 f . ,2 



— / |u(t,j;)|^ dxdr < liminf — / 



(6.1) / \Du\''<C[^) / iDitr VO<p<r<i? 

and similarly, all weak solutions u of^Cu — in Q^{Xq) satisfy 

(6.2) / \Duf<c(^Y^^ [ \Duf yO<p<r<R. 
JqUXo) ^r/ Jq+(Xo) 

Proof. First, we consider the case r = 1 and Xq = 0. Assume that m is a weak 
solution of Cu = in Qi{0). Let k be the smallest integer strictly larger than n/2. 
Since the coefficients A""^ of the operator C are independent of x, we obtain by 
iterative applications of the energy inequalities 

(6.3) \\u{t, ■)\\^.,.(B,M0)) ^ ^ MLHQ-m £ (-1/4, 0), 

where C — C(n, A, A). Moreover, by the Sobolev inequality (see e.g., [8i §5.6.3]) 

(6.4) \\n{t,-)\\^^ 

where C = C{n, A, A, -81/2(0)) = C{n, A, A). Then, we obtain from ([131) and (jO)) 
|M|'<Cp"+2||w||2 Vpe (0,1/2), 

'p(O) ' 



32 



S. CHO, H. DONG, AND S. KIM 



and thus by replacing C by max(2"+^, C) if necessary, we have 

(6.5) / < Vpe(0,l). 

Now, we consider general case. Assume that tt is a weak solution of Cu = in 
Q~{Xq). Then, u{t,x) := u{{t — tQ)/r^,{x — xo)/r) satisfies Cu = in Qi{0), 
where the coefficients £ are still independent of x and satisfies (|2.3|) and (|2.4p : 
more precisely, they are given by A°''^{t) = A°'^{{t — to)/r'^). Therefore, by ()6.5p 
and the change of variables, we obtain 

r-"--' I < C(p/r)"+2 \MIhq,(x.)) Vp G (0, 1). 

The above estimate is equivalent to 

(6.6) [ \u\'<c(^Y^' f \u\' Vpe(0,r). 

JQp{Xo) JQriXa) 

Therefore, (|6.1[) follows from (|6.6p and the observation that for all k = 1, . . . ,rt, 
DkU is also a weak solution of Cu = in (3~(Xo). □ 
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